General Relativistic versus Newtonian: a universality in radiation hydrodynamics 
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We compare Newtonian and general relativistic descriptions of the stationary accretion of self- 
gravitating fluids onto compact bodies. Spherical symmetry and thin gas approximation are as- 
sumed. Luminosity depends, amongst other factors, on the temperature and the contribution of gas 
to the total mass, in both - general relativistic (Lgr) and Newtonian (I/jv) - models. We discover 
a remarkable universal behaviour for transonic flows: the ratio of respective luminosities Lgr/Ln 
is independent of the fractional mass of the gas and depends on asymptotic temperature. It is close 
to 1 in the regime of low asymptotic temperatures and can grow by one order of magnitude for high 
temperatures. These conclusions are valid for a wide range of polytropic equations of state. 

PACS numbers: 



I. INTRODUCTION 

The Newtonian description of accretion has many ad- 
vantages over the general relativistic modelling; it is 
much simpler conceptually, analytically and numerically. 
Therefore it is of interest to specify limits of the validity 
of Newtonian models. It is already known,that they can 
fail in the accretion of transonic test fluids Q-Q. It is not 
suprising that the Newtonian description is not suitable 
for a class of transonic flows in radiation hydrodynam- 
ics. But it is suprising that there emerges a universality 
unknown in the existing literature. We prove that the ra- 
tio Lgr/Ln of general relativistic Lgr and Newtonian 
Ln luminosities is sensitive only to the asymptotic tem- 
perature, assuming the thin gas approximation and some 
additional conditions. 

Models considered here constitute extensions - New- 
tonian and general relativistic - of the classical Bondi 
model [3| to radiation hydrodynamics. They have been 
formulated by several authors (incomplete list includes 
[3 ~ [HI; see also references therein). We make a num- 
ber of assumptions — spherical symmetry, a polytropic 
equation of state and the thin gas approximation in the 
transport equation (.11]. It is assumed that the accretion 
is quasi-stationary. 

There are two kinds of possible general relativistic ef- 
fects. The first is related to the backreaction - and that 
includes self-gravity and dependence on the fractional 
mass content (defined later as 1 — y) of the accreting 
gas - and the other is related to the asymptotic speed of 
sound (we shall occasionally use the term "asymp- 
totic temperature"; both quantities are proportional). 
The same boundary conditions are assumed in both mod- 
els. Each of the two quantities, Lgr and Ln, separately 
depends on y and a^, but their ratio Lqr/Ln depends 
only on the asymptotic temperature. This universality 
means that the enhancement of Lqr/Ln can be found 
by solving the accretion of test fluids in hydrodynamics 
without radiation. This is a much simpler - algebraic - 
problem than the original one. These facts are valid for 



polytropic equations of state p — Kp^ with 1 < T < 5/3. 

The order of the rest of this paper is as follows. Section 
II discusses notation and equations. Section III shows 
how the validity of the thin gas approximation constrains 
the choice of boundary data. We describe the form of 
boundary data for the accretion problem. Section IV 
discusses the notion of sonic points and transonic flows 
in radiation hydrodynamics. In the first part of Section 
V we formulate a low-radiation condition and prove the 
universality of Lgr/Ln- The luminosity of hot accreting 
gases in a general relativistic model can be much higher 
than that given by the related Newtonian counterpart. In 
the second part of Section V we find a sector of luminous 
transonic flows which again gives a universal ratio, but 
now the asymptotic temperature is low and Lgr/Ln ~ 
1. Section VI describes the numerical scheme that is 
used in this work, gives the equation of state and details 
numerical values for most of boundary data. Section VII 
reviews numerical results. The last Section summarizes 
the main conclusions. 



II. FORMALISM AND EQUATIONS 
A. General relativistic accretion. 



The metric 

ds 2 = -N 2 dt 2 + Mr 2 + R 2 (d9 2 + sm 2 {6)d(/) 2 ) (1) 

uses comoving coordinates t, r, < 8 < tt, < <fi < 2tt: 
time, coordinate radius and two angle variables, respec- 
tively. R denotes the areal radius and N is the lapse. 
The radial velocity of gas is given by U — jt^tt- We 



assume relativistic units with G 



1. 



The energy-momentum tensor reads — T E U + T E , 
where the baryonic part is given by T^ v = (p + p) U^U V + 
P9fj,u with the time-like normalized four- velocity U^, 
U^U^ = —1. The radiation part T E V possesses only 

-p E = -TI E and 



four non-zero components, T® E 



2 



T E = T E r . A comoving observer would measure local 
mass densities, the material density p = T B ^ v UnU v and 
the radiation density p E , respectively. The baryonic cur- 
rent reads = PqU^, where po is the baryonic mass 
density. Its conservation is expressed by the equation 



V^" = 0. 



(2) 



Let be the unit normal to a coordinate sphere lying in 
the hypersurface t — const and let k be the related mean 
curvature scalar, k — -^V.jn 1 = -^d r R. The quantity 

j = U^n v NTl lE /VE = NT° E /Vd is interpreted as the 
comoving radiation flux density. We assume the poly- 
tropic equation of state p = Kpy, with constants K and 
r. The internal energy density h and the rest and bary- 
onic mass densities are related by p = po + h, where 
h = p/(T-l). 

There are four conservation equations that origi- 
nate from the contracted Bianchi identities, V^T^ 1 - 8 = 
— V^T^ — F v (here v = 0, r). The radiation force den- 
sity F v describes the interaction between baryons and ra- 
diation. This formulation of general relativistic radiation 
hydrodynamics agrees with that of Park Q, Miller and 
Rezzola Q and (on the fixed, Schwarzschildean, back- 
ground) Thorne et. al @. 

One can find the mean curvature k from the Einstein 
constraint equations G M o = SttT^q ([l|, [l2} ) 



k 



1 - 



2m(R) 
R 



(3) 



where m(R) is the quasilocal mass, 



m{R) = M - 4tt / drr 2 [ p + p E + ) . (4) 



The integration in Q extends from R to the outer bound- 
ary i?oo of the ball of gas. Its external boundary is con- 
nected to the Schwarzschild vacuum spacetime by a tran- 
sient zone of a negligible mass. Thus the asymptotic mass 
M is approximately equal to m (i?oo)- 

In the polar gauge foliation one has a new time 
ts(t,r) with d ts = d t — NUOr. The quantity 

AirNkR 2 (j (l + (f ) 2 ) + 2Up E /k) is the radiation flux 

measured by an observer located at R in coordinates 
(is, R). One can show that 



d ts m(R) = -AirNkR 2 

A-kNUR 2 (p + p) . 




(5) 



The mass contained in the annulus (i?,i?oo) changes if 
the fluxes on the right hand side, one directed outward 
and the other inward, do not cancel. The local baryonic 
flux reads M = —AnU R 2 po, it is not constant {OrM 7^ 0) 
and its boundary value reads M^. 

The accretion process is said to be stationary (or quasi- 
stationary) if all relevant physically observables, that are 



measured at a fixed areal radius R, remain approximately 
constant during time intervals much smaller than the 
runaway instability time scale T = MjM^. That means 
that d ts X = (d t - NUd R )X = for X = p , p, j, U. . .. 

The above assumptions imply that in the thin gas 
approximation Fq = and the radiation force de nsity 
has only one non-zero component F r = nkVdpoj [10| . 
Baryons and radiation interact through the elastic Thom- 
son scattering, n is a material constant, in standard units 
k = aj (m p c) and c, a and m p are respectively the speed 
of light, the Thomson cross section and the proton mass. 

The full system of equations in a form suitable for nu- 
merics has been obtained in 1101. It consists of: 



i) the total energy conservation 



MN- 



r - 1 



r - 1 - a 2 



2MNy- = A11R jNk [ 1 + — ] + C; 



IP 



(6) 

the constant C is the asymptotic energy flux in- 
flowing through the sphere of a radius Roc- 

ii) The local radiation energy conservation 



1 2m ( R ) \ N d I c>2 E\ 

~R / "R 2 "dR > 



dU 



-Kk 2 Njp + 2N (Up E - kj) — + 
2k (jU - kp E ) ^ + 8nNR (j 2 - jp E j ) . (7) 



iii) The relativistic Euler equation (below a — , , 
the speed of sound) 



dR 



r - 1 - a 2 

n2-Ul 



1 fmlR) o , ~ / U 



« j 1 - 



r - 1 



iv) The baryonic mass conservation 
U d 



dU 
dR 



2 2U , 4nRj 
ma 



(8) 



(9) 



r - 1 - a 2 dR "' R ' k 
v) The equation for the lapse 

ifU^..,-^). (10) 



dN f r - 1 - a 2 

dR =N ( KJ — 



Equations ^ and (JS]-[ini) give the complete model 
used in numerical calculations. 

The asymptotic data for the accretion must satisfy 
several physical conditions. We assume the inequali- 
ties a 2 ^ >• M/Roo ^> ensuring, as demonstrated 
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by Karkowski et al. @ and Mach et al. [l3j . that the 
assumption of stationary accretion is reasonably well sat- 
isfied. They are probably required by the demand of sta- 
bility (see a discussion in [9| and studies of stability of 
accreting flows in Newtonian hydrodynamics 13]). In 
the asymptotic region joo rs p^ and the total luminosity 
is well approximated by Lq = 4tt-Rj Oi j 0O . The total lumi- 
nosity is related to the asymptotic accretion rate by 
10] ' 



Lo = aM^|l-^/l-^l 1 M 00 . (11) 



k(R ) 



Ro 



Here Ro is the size of the compact core and the quantity 
a can be interpreted as a binding energy per unit mass. 



B. Newtonian approximation 

The notation is as in the preceding part of this sec- 
tion. The mass accretion flux is now R-independent, in 
contrast to the general relativistic case, <9rM = and 
the baryonic mass density po coincides with p. <j)(R) is 
the Newtonian gravitational potential, 



<f>(R) = -^p--4<n J*~ rp(r)dr; 



(12) 



M(R) — M — Att J r °° r 2 p(r)dr is the mass contained 
within the sphere R. 

The Newtonian model can be described by two basic 
equations @: 

i) the energy conservation equation 



L - L(R) 



M 



r - i 



u 2 

2 



(13) 



ii) The luminosity equation 



L = Lq exp 



-kM 
A-kR 



Lq exp 



— LqRq 

LeR 



(14) 



Notice that the luminosity has the same form as in the 
case of test fluids [IJ . Here we introduced the Eddington 
luminosity L E = A-kM/k while R = GM/\4>(Rq)\ is a 
kind of modified size measure of the compact body. In the 
case of test fluids _R = Ro- We assume Lq = \cj) (Ro) \M; 
notice, however, that for small a this relation (with a = 
\<fi(Ro)\) appears as the Newtonian limit of Eq. (|lll) . 



III. THIN GAS APPROXIMATION AND 
BOUNDARY DATA 

The thin gas approximation demands that the optical 
thicknes [1JJ of the cloud is smaller than one, i.e. 



n(r)adr < 1, 



(15) 



Ho 



where n is the baryonic number density. Notice that 
n = if we assume the monoatomic hydrogenic gas. 
Assuming that n decreases to the asymptotic value rioo, 
we arrive at 1 > RooTi^a. Thus the rough condition 
for the validity of the thin gas approximation is that the 
radiation free path I = l/tn^a) is not shorter than the 
size of the cloud, I > i?oo- This implies po < R P a and 
(taking into account that po « p) estimates the mass of 
gas, M g < j^R 2 ^. Denote the solar mass by M Q and 
define 10 s = ^ 

with the thin gas approximation 



^gf. One obtains an estimate consistent 



Mi 

M 



< 1(T 21 x 10 2s x 



M_ 



TIT- (16) 



We choose s and M that give the right hand side of (jT6|) 
of the order of unity. In such a case a significant part of 
the total mass M would be contributed by the gas itself. 
That could allow for the strong impact of backreaction 
and selfgravitation onto accretion. It is clear that there 
is a scaling freedom - one can trade the size (represented 
by the exponent s) for the total mass without changing 
the bound in (fT6f . 

The boundary data set is the same for the Newtonian 
and general relativistic models. Thus we specify in both 
cases the same values of asymptotic masses M, masses of 
the core, the binding energy per unit mass a ~ \<j) (Ro) |, 
the asymptotic speed of sound and the size R^ . The 
total luminosity Lq is not a free data, but it results from 
equations. We assume identical equations of state in the 
two models. 



IV. SONIC POINTS AND LUMINOSITY 

We shall study transonic flows. For these flows there 
exists a radius i?* such that a* = |J7* |; the speed of sound 
is equal to the length of the spatial part of the velocity 
vector. Henceforth all quantities denoted by asterisk will 
refer to a sonic point. 

It is clear from the inspection of equations that the 
regularity of solutions demands a particular relation for 
the fraction m*/i?*; here m* is the mass within the sonic 
sphere. In the Newtonian model the three characteristics, 
a*Ar, U*n, and / R*n are related as below Q 



l *N 



u. 



Til* jy 



2i? a 



N 



1 - 



N 



III :.\ ( L* N M 



2R, 



N 



Letti 



*N 



(17) 
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In the last equation appears the Eddington luminosity 
Le • It is clear that the necessary condition for the critical 
Newtonian flow - that is, possessing a sonic point - reads 



< 1. 



(18) 



Define x = Lq/ Le and y = m*/M. Since < Lq, the 
inequality x < y becomes the necessary condition for a 
sonic point. In the general relativistic model, at the sonic 
point a 2 = xj-; the denominator of the right hand side 
of Eq. ([5]) vanishes and that implies the vanishing of the 
numerator. One obtains 



1 



k 2 R 

Kj 1 



R- 2U 



47rir p E 



u 

P + J T 



r - 1 



(19) 



Let us remark, that in the Newtonian limit a 2 <C 1 and 
^ r Igr(p*e+P*+Jj:) < l^ff- Therefore, in this 
limit Eq. (fT9]) coincides with Eq. (jTTJ) . It is obvious that 
radiation pushes the sonic point inward; if the size of a 
compact object is bigger than the value of i?* predicted 
by (fT5|) and ([T^| . then the flow becomes subsonic. 



V. UNIVERSALITY IN Lgr/L? 



We assume in this Section the polytropic equation of 
state p = KvF with T < 5/3 — e, for some small e > 0. 
This restriction is due to the peculiar character of the 
equation of state corresponding to T = 5/3 The con- 
stancy of Lgr/Lm is valid for all F's, although the spe- 
cific value of this ratio depends on the equation of state. 



A. Low luminosities 

The natural reference quantity for radiating systems 
is the Eddington luminosity Le- It can be roughly de- 
scribed as the luminosity at which the infall of gas is 
prevented. Thus one might define weakly radiating sys- 
tems as radiating with a luminosity Lq (herein Lq = Lqr 
or Lq = Ljv) that is much smaller than the Eddington lu- 
minosity, Lq <C Le- We will adopt a different definition, 
for reasons that will become clear. 

The (XY) condition. We will say that an accretion 
system satisfies the (XY) condition if a; <C y. 

Notice the trivial fact that y < 1. If (XY) holds, that 
is x <C y, then obviously Lq <C Le- Thus the (XY) 
assumption is stronger than just the statement Lq <C 
Le- Another interesting fact is that (XY) guarantees 
that the characteristics of the sonic point are essentially 
unchanged by the radiation - see Eqs. (fT7|) and (|T9| . 
The luminosity is the product of a by the (asymptotic) 
mass accretion rate, and since the mass accretion rate 
can be formulated completely in terms of the sonic point 



parameters, it becomes luminosity independent if x <C y. 
The general relativistic mass accretion rate Mqr within 
the steadily accreting fluid can be expressed as below (see 
Eq. (6.1) in [l|) 



M, 



GR 



p2 

2 *G R 

m7l *GRPooGR 2 



l *GR. 



l *GR 



*GR 



1 + ?>al GR 



The corresponding Newtonian expression reads 



p2 

TiY 2 ri *N 

M N = Trm* N pcoN—5— 

m *N 



(5-3T) 

2(r-i) 



(20) 



(21) 



Eq. (|2ip has been derived by Kinasiewicz in [14j, but 
it follows also from (|20|) in the limit of small sound 
speeds, ci*gr << 1- The way of writing these two ex- 
pressions is not accidental. It has been shown in [l5[ 



that characteristics of the sonic point 



L *GR 



and 



- do not depend on the fraction of mass carried by the 
gas. These quantities are dictated just by the asymp- 
totic speed of sound aoo in a test fluid model. An anal- 
ogous result holds in the Newtonian model, as shown in 
[14l ]. Therefore Lqr/Ln — Mgr./ Mn is equal to the ra- 
where the coefficient F depends on 



r (and thus on the equation of state) and on the sonic 



point parameters a*N, cl^gr, i GR and 



Therefore 



the coefficient F is independent of the mass fraction y. 
Now the masses are approximately equal, m*GR ~ to*jv; 
this is because the masses within the sonic point are well 
approximated by the masses of the cores, and the lat- 
ter are equal by definition. The equality of masses of 
the cores in both models is one of our boundary condi- 
tions. The asymptotic gas densities PooGR and P™n are 
approximately equal to (M — tti^gr) /V ([3, EH); in 
order to show that one should invoke assumptions con- 



cerning boundary conditions U'^ -M- <C 



The 



calculation is long but straightforward. Thus, we finally 
obtain Lqr/Ln — F; the ratio of luminosities is inde- 
pendent of the fraction of mass carried by the gas, in 
the regime of low luminosities. This means that the ap- 
propriate information on the ratio of the relativistic and 
Newtonian luminosities, Lgr/Ln, can be obtained just 
by the analysis of accreting systems with test gas (and 
for these see, for instance, results in [3 and 0|). This 
is despite the fact that actual values of both luminosi- 
ties taken separately depend on the contribution of the 
gas to total mass. We already know that in accretion 
without radiation the mass accretion rates are maximal 
when m* = 2M/3 and they tend to zero at both ends: 
i) to* — ¥ M (when the density poc tends to zero) and ii) 
to*/M — > (when the mass of the core is negligible in 
comparison to the mass of the fluid) 15]. That implies, 
for weakly radiating systems, that luminosities behave 
in a similar way. But still their ratio is constant and 
independent of the parameter y. 
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B. High luminosities. 

The argument of the former subsection does not apply 
to accreting luminous systems, when the total luminos- 
ity L is close to the Eddington limit or more generally - 
when the (XY) condition is broken. That this can hap- 
pen, is easily illustrated by the Newtonian model. In this 
model one obtains an equation that relates luminosity x 
and the mass content y (see Eq. (21) in in units 
adapted to the convention of this paper): 



5 - 3r 



(5-3r)/(2(r-i)) 



Here Xoo is a constant. The argument that was used 
above relied on the fact that the right hand side of (|22|) 
does not depend on x if x <C y. But if x is relatively large, 
then M becomes x-dependent, and (j2"2l yields a relation 
x = x{y). A similar reasoning can also be applied to the 
general relativistic model; again M depends on a; if a; is 
large. In conclusion, for luminous systems the ratio of 
Lgr/Ln can become x and y-dependent. 

Luminous systems are characterized by small values of 
the asymptotic speed of sound, doo <C 1. We restrict our 
attention to systems that satisfy the following 

X(l-Y) condition. We will say that an accretion 
system satisfies the X(l — Y) condition if x ^> 1 — y and 
x < y/2. 

Since y > x, the above implies y > 2/3. Thus X(l- 
Y) selects a subclass of luminous accretion systems with 
moderate contribution of the gas to total mass. Lumi- 
nous test fluids belong to this category. 

Define an auxiliary quantity 



r P 



L = —2MN— + 4nR 2 jNk 1 + — 
Po ' 



IP 



(23) 



L represents local luminosity as measured by an observer 
stationary at R It follows from Eq. © that at the 
boundary of the accretion cloud L(R oc ) = Lq^; L(R qo ) 
is the total luminosity. By differentiating Eq. ([6]) and 
employing Eq. (1101) one can easily derive the following 
equation 



We have j^M 



-16-K 2 R 3 ^-j It is convenient to 

replace M by = M + k J^°° drjr 3 ^-; the new object 

is constant and coincides with M at the boundary i?oo- 
Then, 



dR 



16n s jN{-R 3p0 



kT-1 



drjr 3 ^ 



Now, the assumptions <C 1 and x < y/2 imply that in 
the region extending from the sonic point to Roa the infall 
velocity U is small and the position of the sonic point i?* 
is large (i?* 3> M). Therefore the metric function k and 
the lapse N are close to unity. The argument bases on 
the approximate validity of Newtonian relations (JTTJ) at 
general relativistic sonic points. The second term on the 
right hand side of (j2"5]l . 



SL = 16n 2 jN 



-R 



3 PO 



F — 1 



fc r — i 



drjr 3 ^- 



is bounded from above by the first term on the right 



( 22 ) hand of (USD, SL < {M^NKj] 



Indeed, the first 
term on the right hand side of (|25|) can be written as 
■4ttMLo jj u ^ fa a ^ - 1S e q ua i to j 47rAIx which in turn is 

larger than jAitMx. The first term of SL is bounded by 
l2n\j\M g = 12irM(l-y). The second term of SL is much 
smaller than the first one if sonic radii are much larger 
than M. Therefore the conclusion follows. 

Thence in the annular region (i?H.,i?oo) the function L 
satisfies with good accuracy the differential equation 



d f »v L 
— L m AiooK -. 

dR 4ttR 2 



(27) 



(25) 



which is solved by L = Lq exp ^— k j^j^j ■ Thus we obtain 
the same form of a solution as in the Newtonian case. 
Furthermore, one can approximate Eq. ([5]) by a suitable 
Newtonian model in the region (R^^R^). Indeed, it is 
easy to show that M w Moo- Expanding the lapse N - 
keeping only the first order terms in a 2 , M(R)/R and U 2 
- we arrive at the Newtonian equation (|13[) . That means 
that under adopted boundary conditions the total general 
relativistic model is well approximated by the Newtonian 
model in the annulus (i?*,i? 00 ), and thus by Eq. (|2"2"|) . 
But Eq. (|2"2"j) has a unique solution, assuming x < y < 1 . 
Therefore the Newtonian limit of the general relativistic 
model and the Newtonian solution do coincide and the 
ratio Lgr/Ln is not only constant, but it is equal to 1. 



VI. NUMERICS. 

We compare two accreting systems, a Newtonian one 
and its general relativistic counterpart, that have inden- 
tical sizes, the same asymptotic masses and identical 
masses of compact cores, equal asymptotic temperature 
and the same binding energy. Thus, it is legitimate to 
say that the boundary data are ultimately: the asymp- 
totic mass, the mass of the core, the binding energy per 
unit mass a = \4>(Rq) |, the asymptotic speed of sound 
doo and the size of the system R^, . The total luminosity 
Lo is not part of these data but is the sought result of 
the two models. 

It appears convenient in numerical calculations to spec- 
ify temporarily pooo and Lq instead of the mass of the 
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core. Conceptually the computational technique is the 
same in the two models. For a given p 0oc one randomly 
chooses Lq (equivalently one could choose an accretion 
rate, due to relation M — Lq/cx). This choice completely 
specifies L(r) in the Newtonian model - see formula (fl4l) . 
Asymptotic radiation data for the general relativistic sys- 
tem in turn are given by jx, = = Lq/ (47ri2^ ), and 
the mass accretion rate = Lo/a. During the numer- 
ical integration one gets a subsonic solution (if the chosen 
Lq is smaller than a critical luminosity) or finds no solu- 
tion at all (if Lq is greater than a critical value). Using 
the bisection method one finds this critical luminosity 
for which the gas flow becomes transonic. The mass of 
the core results from computations. Notice, that for a 
given pooo masses of the core usually differ in the Newto- 
nian and the general relativistic models. The difference is 
particularly noticeable for high asymptotic sound speeds 
Oqc. One should change the value of pooo and repeat the 
procedure until finally the masses of both cores are the 
same for both critical flows. 

In this way one obtains a boundary of the solution 
set (in the plane Lq - M core ) that consists exclusively of 
transonic solutions, if the mean free path of photons is 
larger than the size of the system R^. 

From a mathematical point of view we have a system 
of ordinary first order differential equations. The general 
relativistic problem also includes the integro-algebraic 
constraint Eq. Numerical calculations start from 

the values adopted at the outer boundary Roo and con- 
tinue inward until the equality a = 1 — ffiffl *J 1 — H2lL5i 
(in the GR case) is met at some R; this value of the areal 
radius is denoted as Rq and interpreted as the radius of 
the compact core of the accreting system. For the New- 
tonian model the calculation continues until the gravita- 
tional potential (f> becomes equal to —a. The numerical 
integration employs the 8th order Runge-Kutta method 
fl6j . The main numerical difficulty is encountered in the 
vicinity of the sonic point. In the general relativistic case 
the denominator and the numerator of Eq. (|T9")) vanish 
for a 2 = In numerical computation, the division by 
very small numbers may cause errors and lead to unphysi- 
cal solutions, therefore a special regularization technique 
had to be implemented. We omit further discussion of 
related technicalities, but let us mention that because of 
this difficulty with the sonic point there appear small nu- 
merical errors for M core » 1 and M core ss 0.1 (see Fig. 
3). In the Newtonian model one has to deal with the 
same problem. 

We choose specific numerical data, but since the ac- 
creting system possesses a simple scaling property - as 
discussed in one of preceding sections - one can extend 
the validity of all conclusions to a large family of systems 
with appropriately scaled masses M and sizes i?oo. 

We assume standard gravitational units G = c = 1, the 
polytropic index T = 3/2, the size Roo = 10 s x M and 
the mass M = 1O 6 M , where Mq is the Solar mass. In 
the scaling M = 1 one gets k = 3.62 58 x 10 22 (M©/M), 



that is k — 3.6258 x 10 16 . The Eddington luminosity 
reads L E = 3.4658 x lCT 22 Af/M Q = 3.4658 x 1CT 16 . 
These data are arranged to ensure the validity of the thin 
gas approximation. The optical thickness (1151) is always 
smaller than 1. 



l 




o I 1 , 1 1 — I 

5 -W 9 1 -1(T 8 2 -KT 8 3 -10 1 

Luminosity Lq/Le 



FIG. 1: Luminosity of general relativistic and Newtonian 
models, a = 0.9 and = 10 _1 . The abscissa shows the 
luminosity in terms of the Eddington luminosity Le and the 
ordinate shows the mass of the compact core. 



VII. RESULTS 

Figures[lHl]show accreting solutions on the luminosity- 
(mass of the central core) diagram for a = 0.9. The 
squared speed of sound is a 2 ^ = 10 _1 and a 2 ^ = 10~ 6 , 
respectively. The two figures show transonic solution sets 
for the Newtonian and general relativistic models; they 
are depicted by dashed and solid lines, respectively. Com- 
paring these figures, we notice that the brightness of a 
system increases sharply with the decrease of a 2 ^ . In the 
case illustrated in the second figure, maximal luminosities 
go up to one quarter of the Eddington luminosity. In the 
test gas limit, the interaction between gas and radiation 
is negligible and the gas accretion can be approximated 
by the purely hydrodynamic description. Such a case 
was already analyzed in [l|, with the same conclusion as 
suggested by the comparison of Figs 1 and 2: the larger 
the asymptotic speed of sound, the larger the gap be- 
tween the general relativistic and the Newtonian predic- 
tions. The general relativistic model gives significantly 
larger accretion rates for high asymptotic temperatures. 
These figures clearly demonstrate that luminosities de- 
pend on the fraction of mass deposited in the gas and 
become maximal when this fraction is not bigger than 
1/3. Again, this aspect of the description of the regime 
of weakly radiating sources agrees with the purely hy- 
drodynamic study of [l5| . The position of this maximum 
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FIG. 2: Luminosity of general relativistic and Newtonian 
models, a — 0.9 and a^, = 10 -6 . The abscissa and ordi- 
nate are as in Fig. 1. 



depends weakly on the relative luminosity L/Le and it 
shifts from y = 2/3 in Fig. 1 towards y = 0.75 in Fig. 2. 
It is clear that this effect is due to the influence of the 
radiation; the higher the luminosity, the larger the mass 
of the core at the maximum. 

Fig. 3 reveals a feature of spherical accretion that 
confirms the analytic proof (made in one of preceding 
sections), that Lgr/Ln should be constant, at least 
for small luminosities. While each individual quantity 
Lqr, Ln depends on the contribution of the gas to the to- 
tal mass, their ratio is roughly constant at a given asymp- 
totic temperature. In the six sets of transonic flows (lines 
2-7 on Fig. 4) the ratio Lgr/Ln is independent of the 
mass of accreting gas. We would like to call the reader's 
attention to line 2, where the maximal value of x = 1/4 
is achieved at y = 0.75 (see Fig. 2). Thus the maxi- 
mal value of x/y « 1/3, x ~ 1 — y and still the fraction 
Lgr/Ln is constant. This agrees with the analytic re- 
sult shown in the second part of Section V, although the 
proof of this requires x ^> 1 — y. That suggests that ana- 
lytic results can be proven under less stringent conditions 
than stated in Sec. V. 

Line 1 in Fig. 3 and Fig 4. display data where the 
backreaction effect causes Lgr/Ln to vary (and in par- 
ticular, Lgr/Ln can be made significantly smaller than 
1). Notice, however, that in the general relativistic model 
the flows cease to be transonic for y < 0.94. In contrast, 
they are always transonic in the Newtonian model. There 
is a small segment just below y = 1, where the X(l-Y) 
condition is met and the ratio of luminosities equals 1 . 



VIII. CONCLUSIONS 

There are interesting universal properties hidden in 
jeneralizations of the classical Bondi accretion model. It 



FIG. 3: Binding energy a — 0.9. The values of Lgr/Ln 
are shown on the abscissa. The mass fraction y is put on the 
ordinate. Asymptotic squared speeds of sound are: 10 -7 (line 
no 1); 10~ 6 , 10~ 5 , 10 -4 (the three close lines are grouped as 
line 2); 10" 3 (line 3); 10" 2 (line 4),2.5xl0" 2 (line 5),5xl0~ 2 
(line 6), 10" 1 (line 7). 
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FIG. 4: Binding energy a = 0.9. The luminosities are shown 
on the abscissa while the mass fraction y (the core mass) is 



put on the ordinate. Here a a 



10" 



is already known that when radiation is absent, transonic 
flows (that maximize mass accretion rates) correspond to 
the case when y = m*/M = 2/3, irrespective of the equa- 
tion of state and the asymptotic speed of sound , [3] . 
The mass of the core is about 2/3 of the total mass of 
an accreting system. This paper deals with radiating ac- 
cretion flows. We compare luminosities corresponding to 
transonic solutions of the general relativistic and New- 
tonian accretion models, assuming the same polytropic 
equation of state and identical boundary data - asymp- 
totic speed of sound aoo, size i?oo, total (asymptotic) 
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mass and fraction 1 — y of the total mass contributed by 
gas. We focus our attention on the investigation of the 
relation between their relative luminosity (LgrI ' Ln) and 
y. When accreting systems are characterized by low lu- 
minosity and the condition (XY) of Section V holds true 
(that is Lqr <C Le X y and Ln <C Le x y) then the ra- 
tio Lgr/Ln is independent of y and can be significantly 
larger than 1 . We have found an example with the largest 
value of Lgr/Ln exceeding 1.6, but in earlier investiga- 
tion of test fluids with the polytropic index close to 5/3 
the ratio of mass accretion rates Mqr/M n exceeded 10 



0, which suggests that Lgr/Ln can grow by one order 
of magnitude. On the other hand, when the condition 
X(l-Y) of Section V is valid (that is, a transonic flow 
is highly luminous, x = Lq/Le S> 1 — y and x < y/2, 
but the contribution of gas to the mass is small), then 
Lgr/Ln ~ 1. These properties of the ratio Lgr/Ln 
have been derived analytically and confirmed (under less 
stringent conditions) numerically. 
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